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Abstract 

Background: The circadian clock is an endogenous mechanism that coordinates biological processes with daily 
changes in the environment. In plants, circadian rhythms contribute to both agricultural productivity and 
evolutionary fitness. In barley, the photoperiod response regulator and flowering-time gene Ppd-Hl is orthologous 
to the Arabidopsis core-clock gene PRR7. However, relatively little is known about the role of Ppd-Hl and other 
components of the circadian clock in temperate crop species. In this study, we identified barley clock orthologs and 
tested the effects of natural genetic variation at Ppd-Hl on diurnal and circadian expression of clock and output 
genes from the photoperiod-response pathway. 

Results: Barley clock orthologs HvCCAl, HvGI, HvPRRl, HvPRR37 {Ppd-Hl), HvPRR73, HvPRR59 and HvPRR95 showed a 
high level of sequence similarity and conservation of diurnal and circadian expression patterns, when compared to 
Arabidopsis. The natural mutation at Ppd-Hl did not affect diurnal or circadian cycling of barley clock genes. 
However, the Ppd-Hl mutant was found to be arrhythmic under free-running conditions for the photoperiod- 
response genes HvCOl, HvC02, and the MADS-box transcription factor and vernalization responsive gene Vrn-Hl. 

Conclusion: We suggest that the described eudicot clock is largely conserved in the monocot barley. However, 
genetic differentiation within gene families and differences in the function of Ppd-Hl suggest evolutionary 
modification in the angiosperm clock. Our data indicates that natural variation at Ppd-Hl does not affect the 
expression level of clock genes, but controls photoperiodic output genes. Circadian control of Vrn-Hl in barley 
suggests that this vernalization responsive gene is also controlled by the photoperiod-response pathway. Structural 
and functional characterization of the barley circadian clock will set the basis for future studies of the adaptive 
significance of the circadian clock in Triticeoe species. 
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Background 

The circadian clock is an autonomous oscillator that 
produces endogenous biological rhythms with a period 
of about 24 hours. This clock allows organisms to antici- 
pate predicted daily changes in the environment and to 
coordinate developmental and metabolic processes with 
environmental cues, such as light and temperature, that 
cycle with the rotation of the earth [1-5]. Conceptually, a 
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circadian system comprises three basic components: in- 
put pathways that sense light and temperature, a core os- 
cillator that defines the rhythm, and output pathways 
that control various developmental and metabolic pro- 
cesses, resulting in the optimal adaptation to daily chan- 
ging environments. The core oscillator that generates 
circadian rhythms is comprised of autoregulatory inter- 
locking positive/negative feedback loops. In the eudicot 
Arabidopsis, the central loop called the "core oscillator" 
is composed of two partially redundant Myb-like tran- 
scription factors, CIRCADIAN CLOCK ASSOCIATED 1 
(CCA1) [6] and LATE ELONGATED HYPOCOTYL 
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(LHY) [7], and the pseudo response regulator (PRR) 
TIMING OF CAB EXPRESSION 1 (TOC1). The morn- 
ing expressed CCA1 and LHY repress TOC1 by directly 
binding to its promoter, which results in the evening 
accumulation of TOC1 that in turn represses CCA1 and 
LHY expression [8]. These three genes are critical to sus- 
tain rhythms as the ccal Ihy tocl triple mutant was 
found to be arrhythmic [9]. The core oscillator is further 
fine-tuned by a morning-phased loop and an evening- 
phased loop. The morning loop includes members of 
the pseudo response regulated gene family, PRR3, 
PRR5, PRR7 and PRR9 which contain a pseudo receiver 
domain at the N terminus and a CCT (CONSTANS, 
CONSTANS-LIKE, and TOC1) motif at the C terminus 
[10,11]. PRR transcripts start accumulating after dawn 
sequentially in the order of PRR9, PRR7, PRR5, PRR3, 
and PRR1/TOC1 [12], and it has been shown that PRR7 
and PRR9 repress CCA1 and LHY during the day [10]. 
The evening-phased loop is proposed to include 
GIGANTEA (GI), ZEITLUPE (ZTL), TOC1, and an un- 
known factor Y. GI decreases TOC1 protein level 
through stabilization of the ZTL protein [13,14]. The 
decreased TOC1 protein tends to relieve repression of Y, 
increased Y expression in turn activates TOC1 expres- 
sion, so that Y directly and GI indirectly activate TOC1 
expression [15]. 

Output pathways from the oscillator convey circadian 
rhythms to the various physiological and molecular pro- 
cesses, which include many with agronomic importance, 
such as photosynthesis, growth, phytohormone signaling, 
and photoperiodic flowering [16,17]. The circadian clock 
may thus be a key for improving agronomic performance 
and stress adaptation of crops. Indeed, diurnal expres- 
sion analysis of field grown maize showed that -22% of 
all genes in leaf tissue exhibit diurnal expression patterns 
[18]. In addition, a null allele of the rice GI ortholog 
affected diurnal expression of 75% of all tested genes and 
conferred reduced seasonal adaptability in field grown 
rice [19]. These studies highlight the critical role of the 
clock in cereals. Some orthologs of Arabidopsis core- 
clock genes have been identified in the monocot plants 
rice [20] and Lemna [21,22]. The rice genome was 
reported to encode a single ortholog for LHY and CCA1 
and five PRR orthologs designated as OsPRRl ortholo- 
gous to TOC1, OsPRR73/OsPRR37 corresponding to 
AtPRR7 or AtPRR3, and OsPRRS9/OsPRR9S correspond- 
ing to AtPRRS or AtPRR9 [23]. Over-expression of 
OsCCAl or OsPRRl in Arabidopsis modified circadian 
rhythms [23]. Rice orthologs of TOC1 and PRR 7 partially 
complemented the corresponding Arabidopsis tocl and 
PRR7 mutants, which is consistent with the function of 
these proteins being conserved between monocots and 
Arabidopsis [20,24]. A full compendium analysis of the 
monocot clock awaits to be performed. 



The temperate crop barley, which includes the domes- 
ticated form Hordeum vulgare spp. vulgare and the wild 
subspecies Hordeum vulgare spp. spontaneum, is charac- 
terized by high genetic diversity and good adaptation to 
stress prone marginal environments [25,26]. Interest- 
ingly, adaptation in barley is influenced by the photo- 
period response gene Ppd-Hl, also known as HvPRR37, 
which is orthologous to the rice gene OsPRR37 [23] and 
the Arabidopsis clock gene PRR7 [27]. A natural, reces- 
sive mutation in the CCT domain of Ppd-Hl causes 
photoperiod insensitivity and late flowering in cultivated 
spring barley. In contrast, wild and cultivated winter bar- 
ley genotypes harbor the photoperiod responsive Ppd-Hl 
allele, which induces early flowering under long photo- 
periods. Barley genotypes with a photoperiod responsive 
Ppd-Hl allele are characterized by elevated expression of 
Vrn-H3 (HvFTl) homologous to the Arabidopsis gene 
FLOWERING LOCUS T (FT) [27]. In Arabidopsis, FT is 
the mobile florigen hormone that moves as a protein 
from the leaves through the phloem to the shoot apical 
meristem where it induces the switch from vegetative to 
reproductive growth [28,29]. FT expression is triggered 
by the photoperiod response gene CONSTANS (CO) 
[30]. CO protein is degraded in darkness and expression 
of the protein during the day is crucial for induction of 
the floral activator FT and flowering [30]. It was sug- 
gested that the mutation in Ppd-Hl of spring barley 
delayed flowering time by shifting the diurnal expression 
peaks of the barley CO orthologs HvCOl and HvC02 
into the dark phase, so that the proteins are not synthe- 
sized and Vrn-H3 (HvFTl) not expressed [27]. 

Winter barley is vernalization sensitive, exposure to 
prolonged periods of cold during winter are translated 
into an increased competence to flower in spring. 
Vernalization response is controlled by variation at the 
vernalization genes Vrn-Hl and Vrn-H2 and by the 
MADS -box transcription factors HvVrt2, HvBMl, and 
HvBMlO, which are cereal orthologs of SHORT VEGE- 
TATIVE PHASE (SVP) in Arabidopsis [31,32]. In winter 
barley, Vrn-Hl, with similarity to the Arabidopsis meri- 
stem identity genes APETALA1, CAULIFLOWER, and 
FRUITFUL, is only expressed after exposure to cold [33]. 
Insertions or deletions in the first intron of Vrn-Hl in 
spring barley cause up-regulation of the gene independ- 
ently of vernalization [34]. Spring barley is also charac- 
terized by a deletion of the entire Vrn-H2 locus, which 
includes one truncated and two full sequence ZCCT 
(Zinc finger and CCT domain) genes with no clear 
orthologs in Arabidopsis [35]. In photoperiod-sensitive 
winter barley, Vrn-H2 represses Vrn-H3 (HvFTl) to 
counteract the Ppd-Hl dependent long day induction of 
Vrn-H3 prior to winter. Up-regulation of Vrn-Hl during 
vernalization and consequent down-regulation of Vrn- 
H2 transcript levels in the leaf facilitate the up- 
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regulation of Vrn-H3 during long days mediated by Ppd- 
Hl [36]. HvVrt2, HvBMl, and HvBMlO are floral repres- 
sors, which may act downstream of Vrn-Hl and HvFTl 
in barley, and thus, integrate light and temperature 
dependent regulation of flowering [31,37]. However, the 
effects of variation at Ppd-Hl on circadian expression of 
photoperiod and vernalization response genes have not 
yet been analyzed. The natural mutation in the Ppd-Hl 
gene may affect the photoperiod and vernalization path- 
ways either by changing circadian timing of clock genes 
or by direct control of flowering time genes independ- 
ently from its clock function. 

In this study, we analyzed whether orthologs of Arabi- 
dopsis clock genes are structurally and functionally con- 
served in the temperate crop and long-day plant barley. 
For this we 1) identified barley clock orthologs from 
available genomic databases and 2) analyzed their diurnal 
and circadian expression patterns in two barley geno- 
types differing at the photoperiod response gene and 
clock ortholog Ppd-Hl, We showed that barley clock 
orthologs exhibit a high level of sequence similarity and 
conservation of expression profiles as compared to Ara- 
bidopsis and rice clock genes. The natural mutation at 
Ppd-Hl did not affect expression of clock genes, but 
caused arrhythmicity of clock output genes HvCOl, 
HvC02, and Vrn-Hl under constant conditions. Our 
study provides a characterization of the compendium of 
barley clock genes under circadian conditions, and sets 
the basis to explore the effects of the circadian clock on 
performance in temperate crop species. 

Methods 

Plant material and growth conditions 

The spring barley cultivar Scarlett and an introgression 
line S42IL-107 derived by crossing Scarlett with the wild 
barley accession ISR42-8 were used in this study [38,39]. 
Scarlett has a mutation in the CCT domain of Ppd-Hl and 
is late flowering under LD [27]. The introgression line 
S42IL-107 harbors the photoperiod responsive Ppd-Hl al- 
lele introgressed from wild barley and is early flowering 
under LD. Replicate plants of both genotypes were grown 
in soil in a controlled environment growth chamber (Con- 
viron) at 20°C/18°C day/night with a photoperiod of 16-h 
light and 8-h darkness (LD treatment). After two weeks of 
LD treatment, replicate plants per genotype were har- 
vested for a total of 24-h, starting at the transition to lights 
on in the morning (TO). Night samples were collected in 
the dark. After LD sampling, plants were released into 
continuous light at constant 20°C (LL treatment). Collec- 
tion of leaf material was started at T16 and sampling was 
carried out for a total of 48 -h. Under LD and LL, leaf ma- 
terial was sampled every 2-h, while at the end of the day 
and beginning of the night (or subjective nights) sampling 
was performed every hour. Each sample contains the 



second youngest leaves of three pooled plants. Samples 
were immediately frozen in liquid nitrogen and stored at 
-80°C until required. 

Isolation of clock orthologs 

Sequences of HvGI, Ppd-Hl (PRR37), and PRR73 were 
retrieved from the literature [27,40,41]. To identify ortho- 
logs of CCA1/LHY, PRRl PRR5, and PRR9 we used the 
barley database of DFCI Barley Gene Index Project (TIGR 
gene index project at http://compbio.dfci.harvard.edu/tgi/ 
definitions.html). Expressed sequence tag (EST) sequences 
(and Tentative consensus (TC) sequences,) that had 
high sequence similarity with OsCCAl (Os08g0157600), 
OsPRRl (Os02g40510) and TaPRRl (AK333193), OsPRR59 
(Osllg05930) and OsPRR95 (Os09g36220) were recov- 
ered. Primers were designed on the basis of the homolo- 
gous EST sequences and used to amplify from a cDNA 
pool of Golden Promise a full-length cDNA for HvCCAl 
and a partial gene for HvPRRl (Additional file 1: Table SI). 
RNA extraction and cDNA synthesis were performed as 
described in the next paragraph. Amplifications of 2 uL f 
cDNA were conducted in a 20 uL PCR reaction volume 
containing 0.4 uM of each primer, 1 U of TAQ polymerase 
(GO Taq, Promega), 0.08 mM dNTP, 1.5 mM MgCl 2 . 
Amplification conditions were as follows: 98°C for 5 min, 
35 cycles of 98°C (1 min), 56°C (30 sec) and 72°C (3 min), 
followed by an extension step at 72°C (7 min). Amplified 
fragments were cloned using the TOPO TA cloning kit 
(Invitrogen) following manufacturer instructions. For each 
gene, two independent clones were identified and 
sequenced. Sequence analysis, alignments, and in silico 
translations, were performed using programs within the 
Lasergene® 8 suite (Dnastar, Madison WI). 

Sequence analysis 

Orthologous protein sequences for clock genes in Arabi- 
dopsis and monocots were retrieved from TAIR (http:// 
www.arabidopsis.org/), NCBI (http://www.ncbi.nlm.nih. 
gov/) and Phytozome (http://www.phytozome.net/) data- 
bases. If no protein sequences were found, they were 
generated by virtually translating the mRNA sequences 
from the NCBI database (Lasergene® 8 suite, Dnastar, 
Madison WI). Accession numbers are listed in Table 1. 
For each gene family, sequences were aligned by the 
CLUSTALW method using the MEGA5 program [42]. 
The evolutionary distances were computed using the 
Poisson correction method [43]. All positions containing 
gaps and missing data were excluded. The phylogenetic 
analysis was performed using the neighbor- joining 
method within the MEGA5 package [44]. We used 
10,000 trials to obtain bootstrap values. The alignments 
used for the phylogenetic analysis are available as 
Additional file 2: Figure SI, Additional file 3: Figure S2, 
Additional file 4: Figure S3. 
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Table 1 Genes used in this study (* = translation manually 
corrected; ** = contig not covering the full coding 
sequence, first part of the protein is missing) 



Species 


Gene name Gene ID 


Protein ID 


Arabidopsis 
thaliono 


ArPRRU 

/\ir mm / / 

TOO 


ATSG61 380 


AT5G61380 




AtPRR3 


AT5G60100 


AT5G60100 




AtPRR7 


AT5G02810 


AT5G02810 




AtPRR5 


AT5G24470 


AT5G24470 




AtPRR9 


AT2G46790 


AT2G46790 




AtGI 


AT1G22770 


AT1G22770 




AtCCAl 


AT2G46830 


AT2G46830 




AtLHY 




AT1G01060 


Oryza sotiva 


OsPRR 1/ 
TOO 


LOC Os02g40510 


LOC_Os02g40510 




OsPRR37 


LOC_Os07g49460 


LOC_Os07g49460 




OsPRR73 


LOC_Os03g 17570 


LOC_Os03g 17570 




OsPRR59 


LOC_Os11g05930 


LOC_Os11g05930 




OsPRR95 


LOC_Os09g36220 


LOC_Os09g36220 




OsGI 


NM_001 048755 


Os01g01 82600 




OsCCA 1 




NP_001061032 


Sorghum 
bicolor 


ShPRR 1 / 

uur mm // 

TOO 




Sb04g026190 




SbPRR37 


Sb06g014570 


Sb06g014570 




SbPRR73 


Sb01g038820 


Sb01g038820 




SbPRR59 


Sb05g003660 


Sb05g003660 




SbPRR95 


Sb02g030870 


Sb02g030870 




SbGI 


Sb03g003650 


Sb03g003650 




SbCCA 1 


XM 007443R4S 


XP_002443890* 


Zea mays 


7mPRR1/ 
TOO 


mm nm 1 S43S1 


NP_001 147823 




ZmPRR37 


LOC 100280240 


N POO 1146641 




ZmPRR73 


EU9521 16 


ACG24234 




ZmPRR59 


GRMZM2G1 35446 


GRMZM2G 135446 




ZmPRR95 


NM_001 158064 


NP_001 151536 




ZmGI 


BK006299 


DAA06172 




ZmCCAl 


LOC1 001 92868 


NP_001 131529 


Triticum 
oestivum 


ToPRRl/ 
TOO 


AK333193 


in silico 




ToPRR37 


DQ885766 


ABL09477 




ToPRR73 


Contig built from the 

following TC: 

G11 8.1 11D24F01 0720, 

TC37276, TC376302, 

TC377674, 

TC391788,TC392931, 
TC400719,TC418129, 
TC434756 


in silico 




TaPRR59 


Contig built from the 


in silico** 



following TC: BJ298369, 
G608.119D1,TC376876, 
TC393119,TC399029, 



Table 1 Genes used in this study (* = translation manually 
corrected; ** = contig not covering the full coding 
sequence, first part of the protein is missing) (Continued) 



Tritin im 

i 1 1 u v_m i i i 

aestivum 


TnPRR^Q 

I LAI 1 1 / \J s 


TC411244, TaE05039B10R, 
WHE2989, wr1.pk0105.g4 






TaPRR95 


Contig built from the in silico 
following TC: FGAS012265, 
TC373568, TC436795 




TqGI 


AF543844 


AAQ1 1738 




TaCCAl 


HQ222606 


ADW09013* 


Hordeum 
vulgore 


HvPRRl/ 
TOO 


JN603243 






HvPRR37 


AY970701 


AAY42109 




HvPRR73 


AK376549 


BAK07744 




HvPRR59 


AK361360 


BAJ92567 




HvPRR95 


AK252005 


in silico 




HvGI 


AY740524 


AAW66946 




HvCCAl 


JN603242 




Brachypodium 


BdPRRl/ 
TOO 


Bradi3g48880 


Bradi3g48880 




BdPRR37 


Bradi1g16490 


Bradi1g16490 




BdPRR73 


Bradi1g65910 


Bradi1g65910 




BdPRR59 


Bradi4g24970 


Bradi4g24970 




BdPRR95 


Bradi4g36077 


Bradi4g36077 




BdGI 


Bradi2g05230 


Bradi2g05230* 




BdCCAl 


Bradi3g16510 


Bradi3g16510* 



RNA extraction, cDNA synthesis and real time qRT-PCR 

Total RNA was extracted from 100 mg of tissue using 
TRIZOL® reagent (Invitrogen) following manufacturers 
instructions, except for the addition of RNaseH, followed 
by a DNase treatment (final volume 100 uL). First strand 
cDNA synthesis was performed on 4 uL of total RNA 
using 100 U of SuperScriptTM II RT (Invitrogen) and 
500 ng of poly-T primer and following manufacturers 
recommendations (final volume 40 uL). The resulting 
cDNA was diluted 1:4 in nuclease-free water and stored 
in aliquots at -20°C 

Real-Time-quantitative PCRs (qRT-PCR) were per- 
formed on cDNA samples using gene-specific primers 
(Additional file 1: Table SI). Amplifications were per- 
formed using 4 uL of cDNA, 0.5 U of GoTaq DNA poly- 
merase (Promega), 0.2 mM dNTP, 2.5 mM MgCl 2 , 
0.2 uM each primer, and 0.5 uL of EvaGreen (Biotium) 
in a final volume of 10 uL. Reactions were performed in 
a LightCycler480 (Roche) with the following amplifica- 
tion conditions: 95°C for 5 min, 45 cycles of 95°C (10 s), 
60°C (10 s) and 72°C (10 s). Appropriate non-template 
controls were included in each 384-well PCR. Dissociation 
analysis was performed at the end of each run and the 
melting curves for each primer pair showed a single peak 
confirming the specificity of the reaction. The standard 
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curves were prepared from a dilution series of plasmids 
containing the target fragments and subjected to qRT-PCR 
analysis with the respective cDNA samples. Starting 
amounts for each data point were calculated based on the 
titration curve for each target gene and the reference 
(HvActin) gene using the LightCycler 480 Software (Roche; 
version 1.5). Expression values shown in the figures repre- 
sent the average ± standard deviation of 2 technical repli- 
cates of the ratio between the target and reference gene 
values. 

Results 

Phylogenetic analyses of the barley clock genes 

Barley sequences with high similarity to OsCCAl, OsPRRl, 
OsPRR59, and OsPRR95, orthologous to CCA1/LHY, PRR1, 
PRR5 and PRR9 in Arabidopsis were recovered from the 
barley EST databases. Based on the tentative consensus 
sequences (TCs), we obtained a full-length clone for 
HvCCAl (accession number JN603242) and a partial se- 
quence for HvPRRl (accession number JN603243) from 
Golden Promise cDNA. Barley gene sequences of HvGI 
[40], HvPRR37 (Ppd-Hl) [27], HvPRR73 [41] were retrieved 
from the gene banks. In order to examine the structural 
conservation of these potential circadian clock genes in bar- 
ley, candidate clock sequences from barley were compared 
in multiple protein alignments to orthologous sequences 
from Arabidopsis and monocot species, including rice, Bra- 
chypodium, maize, sorghum, and wheat (Additional file 5: 
Table S2). 

In barley, five distinctive PRR sequences could be identi- 
fied. These fell into three major clades, the PRR1/TOC1 
clade, the PRR3/7 clade, and the PRR5/9 clade (Figure la). 
Orthologs of OsPRRl clustered with PRR1/TOC1 from Ara- 
bidopsis and the PRR 1 -like genes could clearly be recog- 
nized as an outgroup of the PRR gene family. Each of the 
two remaining PRR clades showed two subgroups for 
monocot orthologs OsPRR37, OsPRR73 and OsPRR59, 
OsPRR95, respectively. AtPRR3 and AtPRR7 were outgroups 
to the two monocot gene clusters PRR37 and PRR73. 
AtPRRS and AtPRR9, however, grouped with the monocot 
PRR95 genes, while the monocot PRR59 sequences formed 
a separate cluster. Genetic distances were calculated from 
protein sequences for entire orthologous gene sequences or 
for domains within the monocot species and across the 
monocot species and Arabidopsis (Table 2). Variation in 
genetic distances between domains, genes and group of 
genes may indicate differences in evolutionary history and 
in the conservation of gene functions. Within the PRR 
genes, average genetic distances across Arabidopsis and 
monocots were lowest for the PRR1 gene group with a gen- 
etic distance of 0.39, followed by PRR3/7 with 0.51 and 
PRR5/9 with 0.86 (Table 2). The derived PRR amino-acid 
sequences were found to be most similar within their 
pseudo-receiver domain (mean genetic distance 0.32), and 



at the C-terminal end (mean genetic distance 0.22), which 
contains a CCT motif that is commonly also found in the 
CONSTANS (CO) family of flowering regulators [45] 
(Table 2, Additional file 4: Figure S3). 

The identified HvCCAl protein sequence showed a 
high similarity to OsCCAl and the redundant pair 
CCA1/LHY in Arabidopsis, as evidenced by the phylo- 
genetic tree (Figure lb). Sequence distances between 
HvCCAl and OsCCAl were 0.26, while the overall mean 
distance across CCA1 and LHY in Arabidopsis and the 
monocot species was 0.55 (Table 2, Additional file 5: 
Table S2). The CCA1 protein sequences exhibited the 
highest sequence identity of 0.95 at the single MYB 
DNA-binding domain at its N-terminal end (Table 2, 
Additional file 3: Figure S2) [6,7]. 

Multiple sequence alignments using the GI proteins 
from the six monocot plants and Arabidopsis showed 
that the GI protein was characterized by the lowest over- 
all mean sequence distance (0.18, Table 2, Figure lc). 
HvGI was most similar to TaGI and BdGI with sequence 
distances of 0.03 and 0.07, respectively (Additional file 2: 
Figure SI, Additional file 5: Table S2). GI thus appears to 
be evolving slowly. 

From our alignments we concluded that the barley 
sequences HvCCAl and HvPRRl, and the identified clones 
AK376549 (HvPRR73), AK361360 (HvPRRS9) and 
AK252005 (HvPRR95) were the likely orthologs of CCA1/ 
LHY, TOCl/PRRl PRR3/7 and PRRS/9 from Arabidopsis. 
In addition, mean distances calculated for protein 
sequences of the different clock genes across six different 
monocot species revealed that GI showed the highest con- 
servation with a mean distance of 0.10, followed by PRR1 
(0.22), and CCA1 (0.27). The two clades of PRR37/73 and 
PRR59/95, which have evolved after the split of the eudicot 
and monocot lineages, showed the highest mean distances 
with 0.46 and 0.90, respectively. 

Diurnal and circadian expression patterns of putative 
barley clock orthologs 

In order to study the functional conservation of tran- 
script accumulation rhythms in the circadian clock of 
barley, clock orthologs were analyzed for diurnal and cir- 
cadian expression patterns from young leaf tissue. In 
addition, we tested the effects of the natural mutation in 
the CCT domain of Ppd-Hl (PRR37) in the spring barley 
Scarlett on diurnal and circadian expression of clock and 
clock output genes. Plants of both genotypes were 
entrained under long day photoperiods (16-h light/8-h 
dark) for 2 weeks. Subsequently, leaf samples were har- 
vested every two hours under light/dark (LD) for 24 
hours and under continuous light (LL) for 48 hours. 
Additional sampling (every hour) was performed at the 
end of the day and beginning of the night (or subjective 
nights). This provided both a diurnal and free-running 
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Figure 1 Phylogenetic relationships of circadian clock related proteins. Neighbor-joining phylogenetic trees were built from the protein 
alignments of the indicated sequences, a) PRR gene family, b) CCA1/LHY gene family, c) Gl genes. Accession numbers are given in Table 1. 
Bootstrap values were calculated from 10,000 permutations. 



sampling series. The expression of circadian-clock ortho- 
logs in barley seedlings oscillated under LD, and the 
rhythm was sustained under LL conditions. The PRR 
orthologs were expressed in a sequential manner in the 
order of HvPRR73/HvPRR37, HvPRR95/HvPRR59 and 
HvPRRl. Expression of HvPRR3 7 (Ppd-Hl) and HvPRR73 
started rising at T2, HvPRR95 at T4, HvPRR59 at T6, and 
HvPRRl at T8. Expression peaks of HvPRR37 (Ppd-Hl) 
and HvPRR73 were broader than those of the other 
HvPRR genes, so that expression peaks of HvPRR37 (Ppd- 
Hl) and HvPRR73 coincided with that of HvPRR59 and 
HvPRR95 at T8, whhile HvPRRl peaked at T12 (Figure 2). 



HvCCAl expression started rising before subjective dawn 
and showed rhythmic peaks 4 hours after light on under 
LD and 4 hours after subjective dawn under LL 
(Figures 3a, b). HvGl exhibited an expression peak be- 
tween T10-T14 under LD conditions and T12-T14 under 
LL conditions (Figures 3c, d). We did not observe signifi- 
cant differences in the expression patterns, period or amp- 
litude of barley circadian-clock genes between the two 
genotypes differing at the natural mutation in Ppd-Hl. 
The sustained oscillations of barley clock genes under 
free -running conditions showed that these are under circa- 
dian control (Figure 2, Figure 3). The natural mutation in 
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Table 2 Genetic distances for gene families across 
monocots and Arabidopsis for entire sequences or 
sequence domains 



Sequences 


N° of 

sequences 


Domains 


Overall mean 
distances 


Gl a 


7 


InULIcciI IULdllZ.cU piULclll 


0 1 8 


Gl (monocots) b 


6 


Nuclear localized protein 


0.10 


CCA1/LHY 


8 


MYB transcription factor 


0.55 


CCA1 (monocot) 


6 


MYB transcription factor 


0.27 


AtCCAl and AtLHY 


7 


MYB transcription factor 


0.45 


CCA I and AtLLAl 


7 


MYB transcription factor 


0.44 


CCA1/LHY: MYB 
domain 


8 


MYB domain 


0.05 


CCA1/LHY: MYB 
domain (monocot) 


6 


MYB domain 


0.01 


PRRs 


35 


pseudoreceiver-CCT 
domain 


0.87 


PRR1 


7 


pseudoreceiver-CCT 
domain 


0.39 


PRR1 (monocots) 


6 


pseudoreceiver-CCT 
domain 


0.22 


nnno/7 
rKno/ / 


14 


pseudoreceiver-CCT 
domain 


U.J 1 


PRR3/7 (monocots) 


12 


pseudoreceiver-CCT 
domain 


0.46 


PRR5/9 


14 


pseudoreceiver-CCT 
domain 


0.86 


PRR5/9 (monocots) 


12 


pseudoreceiver-CCT 
domain 


0.90 


PRRs: REC domain 


35 


Receiver domain 


0.32 


PRRs:REC domain 
(monocots) 


30 


Receiver domain 


0.31 


PRRs:CCT domain 


35 


CCT domain 


0.22 


PRRs:CCT domain 
(monocots) 


30 


CCT domain 


0.22 



a Arabidopsis and monocots. 
b only monocots. 



the CCT domain of Ppd-Hl in barley did not affect 
diurnal or circadian expression of barley clock gene 
orthologs. 

Diurnal and circadian expression of clock output genes 

In Arabidopsis, important circadian-controlled pro- 
cesses include the photoperiod-dependent control of 
flowering time [46,47]. In order to further characterize 
the barley circadian clock and test for effects of the 
natural mutation at Ppd-Hl on clock output pathways, 
we analyzed diurnal and circadian expression of candi- 
date genes from the photoperiod pathway. We selected 
HvCOl and HvC02, the barley orthologs of the circadian 
controlled photoperiod response gene CONSTANS 
[41,48] and HvFTl, the barley ortholog of the Arabidopsis 
FT [49]. 



In addition, three MADS box transcription factors 
involved in flowering time control in barley were 
selected. These include Vrn-Hl, a flowering activator pu- 
tatively downstream of HvFTl [32], and two repressors 
of flowering, HvVRT2 and HvBMl orthologous to SVP 
(Short Vegetative Phase) in Arabidopsis [31]. Expression 
profiles under LD and LL conditions showed that all six 
genes cycled under diurnal and circadian conditions in at 
least one genotype (Figure 4). Differences in expression 
between the two genotypes differing at Ppd-Hl were 
observed for HvCOl, HvC02, HvFTl, Vrn-Hl, HvVRT2, 
and HvBMl. Under LD conditions, HvCOl and HvC02 
expression cycled in both genotypes, but HvC02 expres- 
sion was lower in Scarlett {Ppd-Hl) than S42-IL107 
(Ppd-Hl) between T16 and T18 (Figures 4a, c). Under 
LL, HvCOl and HvC02 continued to cycle approxi- 
mately once every 24-h in S42-IL107, while expression 
peaks of both genes were two- to threefold lower in 
Scarlett as compared to S42-IL107 (Figures 4b, d). In 
addition, in Scarlett, HvCOl expression peaked at T18 
and T46 and showed a strong delay in oscillation peaks 
under LL. HvC02 only displayed a circadian rhythm in 
S42-IL107 (Ppd-Hl) and showed an arrhythmic and low 
expression under LL in Scarlett, which harbors the 
mutated Ppd-Hl allele. HvFTl expression was below de- 
tection in Scarlett, but in S42-IL107, it oscillated with a 
peak in the afternoon (T13), both under LD and LL 
(Figures 4e, f). Vrn-Hl showed a significantly lower ex- 
pression in Scarlett than in S42-IL107 under LD 
(Figure 4g). Under LL, Vrn-Hl expression only cycled in 
S42-IL107, but not in Scarlett (Figure 4h). HvVRT2 and 
HvBMl cycled under LD and LL conditions in both gen- 
otypes (Figures 4i-l). S42-IL107 (Ppd-Hl) exhibited a 
lower expression of HvVRT2 than Scarlett (Ppd-Hl) 
under LD and LL. All analyzed flowering-time genes 
were thus characterized by rhythmic expression that per- 
sisted under LL conditions. However, the natural muta- 
tion in Ppd-Hl strongly delayed and dampened circadian 
rhythms of HvCOl, HvC02 and Vrn-Hl under LL. Fur- 
thermore, HvFTl expression could not be detected in 
Scarlett with the mutated Ppd-Hl allele. 

In addition, we analyzed diurnal and circadian expres- 
sion of the chlorophyll a/b-binding protein HvCABIII 
[50] and HvGRP7 [51], the barley ortholog of the Glycine 
Rich RNA-binding protein GRP7, also termed cold- 
circadian rhythms-RNA binding (CCR2) [52]. These two 
genes are known clock output genes from the photosyn- 
thetic pathway [53] and slave (non-self-sustaining) oscil- 
lator [51], respectively. We found that both genes 
showed rhythmic expression under LD and LL and no 
significant differences were observed between genotypes 
(Figures 5a-d). These two genes were thus under circa- 
dian control in barley and were not affected by variation 
at Ppd-Hl. 
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Figure 2 Diurnal and circadian expression of PRR genes in the spring cultivar Scarlett ippd-Hl) and the introgression line S42 IL-107 
(Ppd-Hl) under long day and free running conditions. Scarlett (solid line) and S42-IL107 (dashed line) were grown under long day conditions 
(16 hours light - left panel) for two weeks and then released in continuous light (right panel). Transcript accumulation was measured at two-hour 
intervals (one-hour interval at the end of the day and beginning of the night/subjective nights) by qRT-PCR analysis of specific genes and 
normalized to HvActin. Values represent the average of two technical replicates ± standard error. Black and grey bars indicate objective and 
subjective nights. 
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Discussion 

The structure and expression profiles of barley clock 
orthologs are conserved 

Early work in wheat has already demonstrated that 
monocots possess an endogenous circadian clock, Cab-1 
gene expression continued to cycle in wheat plants that 
had been transferred to continuous light or darkness 
[54]. The recent availability of monocot plant-genome 
sequences allows determination of the extent to which 
the genetic clock model developed from Arabidopsis is 
conserved and can thus explain endogenous cycling of 
gene expression in monocot crops. In this work, we 
demonstrated structural similarities between the poten- 
tial barley clock orthologs and rice and Arabidopsis clock 
genes (Figure 1, Additional file 2: Figure SI, Additional 
file 3: Figure S2, Additional file 4: Figure S3, Additional 
file 5: Table S2). PRR genes from six different monocot 
species and Arabidopsis as an eudicot species clearly 
separated into three clades (PRR1/TOC1 clade, PRR3/7 
clade, and PRR5/9 clade) (Figure 1). Each clade con- 
tained genes from both eudicots and monocots, confirm- 
ing that ancient PRR gene(s) diverged into three clades 
before the speciation of monocots and eudicots [55]. 
Within the PRR3/7 and PRR5/9 clade the relationships 
between the Arabidopsis genes and monocot orthologs 
could not be resolved confirming that since the diver- 
gence of monocot and eudicots these clades have 
expanded independently in both lineages as a result of 
genome duplication [55]. Additionally, the high average 
distance of PRR37/73 and PRR59/95 orthologs, as com- 
pared to the other clock genes, suggests functional diver- 
gence and specialization that occurred in the monocot 
lineage. In contrast, CCA1 sequences showed a high 
similarity within the monocots (0.27), as compared to 
the two paralogous genes, CCA1 and LHY in Arabidop- 
sis (0.75). This suggested conservation of CCA1 gene 
function in monocots as compared to stronger diversifi- 
cation between the paralogous genes in Arabidopsis 
(Figure lb, Table 2, Additional file 5: Table S2). Two 
CCA1 paralogs were found in the monocotyledonous 
plant Lemna, suggesting that duplications/deletions of 
clock genes occurred throughout the evolution of eudi- 
cots and monocots [56]. Phylogenetic analysis of GI 
sequences confirmed high amino-acid sequence iden- 
tities within the monocots and between monocots and 
Arabidopsis. High levels of conservation of GI have 
already been demonstrated for seed plants, including 
monocotyledonous plants, such as rice (Oryza sativa) 
[57], Brachypodium [58], wheat [59], barley [40], and the 
eudicotyledons, such as pea [60]. Our data suggest that 
Arabidopsis clock genes are conserved in barley and 
other monocot species. However, differences in gene 
number (CCA1/LHY) and genetic differentiation within 
gene families (PRR gene family) suggest evolutionary 



modification of clock related components between barley 
and Arabidopsis. 

The similarities between the clock amino-acid sequences 
of H. vulgare and Arabidopsis suggest similar molecular 
functions of these proteins. Barley clock genes cycled 
under free-running conditions and their temporal expres- 
sion profiles were similar to those in Arabidopsis. This 
indicated that transcriptional regulatory mechanisms are 
likely conserved between these species. Functional conser- 
vation of some clock orthologs has been demonstrated for 
the monocotyledonous plants rice [21] and Lemna [22]. 
Conservation of expression patterns of barley clock ortho- 
logs under LD and LL support the suggestion that clock 
orthologs function in the circadian clock in barley. 

Natural genetic variation in Ppd-Hl affects photoperiod 
output genes 

Despite structural and functional conservation, diver- 
gence between Arabidopsis and monocot clock genes is 
suggested by evolving functions in the PRR7 lineage in 
monocots. In the monocots, barley and wheat, PRR37 
(Ppd-Hl and Ppd-Al-Dl, respectively) is one of the 
major determinants of photoperiod sensitivity and flow- 
ering time [27,61]. In Arabidopsis, PRR7 contributes to 
the determination of flowering time, although the effects 
are not large and PRR7 may not be major determinant of 
flowering time among natural populations [62]. In Arabi- 
dopsis expression of the PRR1/TOC1 family members is 
under the control of a coordinate circadian rhythm at 
the level of transcription such that the PRRs mRNAs 
start accumulating sequentially after dawn with 2 to 3 
hours intervals in the order of PRR9, PRR7, PRRS, PRR3, 
and PRR1/TOC1 [63]. It was demonstrated that light 
stimuli and possibly phytochromes are crucial to induce 
the early transcription of PRR9 } and this event appears to 
link the intrinsic oscillation of the AtPRRl/TOCl circa- 
dian waves to the entrained external photoperiod [64]. 
Like in Arabidopsis, PRR orthologs in barley were also 
expressed in a sequential manner at approximately 2 
hours intervals from each other (Figure 2). However, 
in contrast to Arabidopsis Ppd-Hl (HvPRR37) and 
HvPRR73, and not the barley orthologs of PRR9, showed 
the earliest increase in expression after dawn and were 
followed by expression of HvPRR95/HvPRR59 and 
HvPRRl. This expression pattern was similar to the 
sequential expression of PRRs in rice with OsPRR73 
(OsPRR37) induced first after dawn, followed by 
OsPRR9S (OsPRRS9) and OsPRRl [23]. In rice, the 
circadian-shape of OsPRR73 was markedly affected 
by the photoperiodic conditions, whereas OsPRRl 
was not affected by changes in photoperiod. Such 
photoperiodicity-dependent changes in the circadian- 
profiles have been reported for certain circadian con- 
trolled Arabidopsis genes, which include the CO gene 
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Figure 3 Diurnal and circadian expression of HvCCAl and HVGI in the spring cultivar Scarlett (ppd-HI) and the introgression line 
S42-IL107 (Ppd-HI) under long day and constant light. Scarlett (solid line) and S42-IL1 07 (dashed line) were grown under long day conditions 
(16 hours light - left panel) for two weeks and then released in continuous light (right panel). Transcript accumulation was measured at two hour 
intervals (one-hour interval at the end of the day and beginning of the night/subjective nights) by qRT-PCR analysis of specific genes and 
normalized to HvActin. Values represent the average of two technical replicates ± standard error. Black and grey bars indicate objective and 
subjective nights. 



that encodes a crucial regulator of flowering time. 
These findings suggest differences in the control of 
PRRs between Arabidopsis and the monocots rice and 
barley. Differential regulation of PRRs in Arabidopsis 
and monocots may be the key to understand why nat- 
ural mutations in Ppd-HI (HvPRR37) and (Hd2) 
OsPRR73 cause variation in flowering time in barley 
and rice [27,65], while mutations of PRR7 do not have 
strong effects on development in Arabidopsis. 

The PRR7 orthologs in barley and wheat may affect 
flowering through changing circadian parameters or 
through a clock independent function of the gene. In order 
to pursue this question, we tested the effects of a natural 
mutation at Ppd-HI on diurnal and circadian expression 
profiles of core clock and output genes. Expression of bar- 
ley clock genes and the clock output genes HvCABIII and 
HvGRP7 did not differ between the two genotypes harbor- 
ing a different Ppd-HI allele under LD or LL (Figures 2, 3, 
5). These results suggest that the mutation in the CCT do- 
main of Ppd-HI does not affect diurnal and circadian cyc- 
ling of clock genes in barley. In contrast, expression 
profiles of barley flowering-time genes showed significant 
differences between both genotypes. The Ppd-HI mutant 
was found to be arrhythmic under free-running conditions 
for the photoperiod-response genes HvCOl, HvC02, and 



the MADS -box transcription factor and vernalization re- 
sponsive gene Vrn-Hl (Figure 4). A wild- type Ppd-HI al- 
lele was thus necessary to maintain circadian oscillations 
of these genes under constant light. Cycling of HvCOl, 
HvC02, and Vrn-Hl was observed in Scarlett {Ppd-HI) 
under LD as compared to LL. This suggested that light/ 
dark cues were necessary for diurnal oscillations in the 
presence of a mutated Ppd-HI allele. Separate effects of 
the mutation in Ppd-HI on circadian expression of flower- 
ing time genes and clock orthologs indicate independent 
functions of this gene in the clock and photoperiod 
pathways. 

PRR9, PRR7, and PRR5 positively regulate flowering 
time in Arabidopsis through the activation of CO expres- 
sion during daytime [66]. Our results show that under 
LL, variation at Ppd-HI also controlled HvCO-lihe gene 
expression in barley. However, under LD, variation at 
Ppd-HI did not show clear effects on diurnal expression 
profiles of HvCOl and HvC02; while the mutation in 
Ppd-HI had a significant effect on HvFTl expression 
levels also under LD. Ppd-HI may thus also affect HvFTl 
expression independently of HvCO-lihe gene transcript ac- 
cumulation under LD. These results are supported by a re- 
cent analysis of a spring barley Bowman and a derived 
introgression line carrying a mutation in HvEl/3, the barley 
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(See figure on previous page.) 

Figure 4 Diurnal and circadian expression of flowering genes in the spring cultivar Scarlett {ppd-Hl) and the introgression line S42 
IL-107 {Ppd-Hl) under long day and free running conditions. Scarlett (solid line) and S42-IL107 (dashed line) were grown under long day 
conditions (16 hours light - left panel) for two weeks and then released in continuous light (right panel). Transcript accumulation was measured at 
two hour intervals (one-hour interval at the end of the day and beginning of the night/subjective nights) by qRT-PCR analysis of specific genes 
and normalized to HvActin. Values represent the average of two technical replicates ± standard error. Black and grey bars indicate objective and 
subjective nights. 



ortholog of Elf3 in Arabidopsis. Elevated levels of HvFTl 
expression and early flowering in the introgression line as 
compared to the recurrent parent Bowman did not correl- 
ate with higher expression of HvCOl [67]. It has already 
been shown that functional variation at Ppd-Hl and Ppd- 
Al-Dl has a major effect on expression of HvFTl and 
TaFT m barley and wheat, respectively [27,61]. Structural 
conservation of PPD1 in wheat and barley (Figure 1) may 
thus reflect also functional conservation. However, in con- 
trast to barley, functional variation at Ppd-Dla in wheat is 
associated with a deletion in the promoter and mis- 
expression of the gene, in contrast to Ppd-Hl in barley 
that harbors a change in the protein coding sequence 
[27,61]. 

The PRR37 ortholog of the short day plant Sorghum 
was recently identified as a repressor of flowering under 



long day conditions [68]. The repressive function of 
SbPRR37 was due to a bimodal expression with peaks in 
the morning and evening which were controlled by light 
and the circadian clock. Loss-of-function mutations in 
SbPRR37 and short day abolished the evening phase ex- 
pression of SbPRR37 and caused significant effects on 
flowering time and expression of SbCO and SbFT, the 
Sorghum orthologs of CO and FT, Variation at both 
genes, Ppd-Hl and SbPRR37 showed significant effects 
on expression of the FT orthologs and flowering in the 
respective inductive day-length. In contrast to the study 
on Sorghum, we show that the floral inducer Ppd-Hl in 
the long-day-grass barley does not exhibit the light 
dependent second expression peak in the evening under 
LD. In addition, the hypomorphic mutation in Ppd-Hl 
caused a reduced and arrhythmic expression of HvCOl/ 



7.0 

c 6.0 

| 5.0 

J 4.0 

5 3.0 
5 



1.0 
0.0 





























ft p 




7.0 

.S 6.0 
o 

^ 5.0 
J 4.0 
§3.0 
I 2.0 
1.0 
0.0 



1 



10 





n 












n 








*. \ 


1 


? / x 






CD O 
-i- OJ 
\— \— 


1 

■nT 00 OJ CO 
C\J OJ CO CO 
— 1 — 1 — 1 — 


1 

o <^ 
h- h- 


■ 

CO OJ CD O ^ 
■^1" LO LD CD CD 
1— 1— 1— 1— 1— 


































1 


1 1 


■ 



Figure 5 Diurnal and circadian expression of clock output genes HvCABII and HvGRP7 in the spring cultivar Scarlett (ppd-Hl) and the 
introgression line S42 IL-107 {Ppd-Hl) under long day and free running conditions. Scarlett (solid line) and S42-IL1 07 (dashed line) were 
grown under long day conditions (16 hours light - left panel) for two weeks and then released in continuous light (right panel). Transcript 
accumulation was measured at two hour intervals (one-hour interval at the end of the day and beginning of the night/subjective nights) by qRT- 
PCR analysis of specific genes and normalized to HvActin. Values represent the average of two technical replicates ± standard error. Black and grey 
bars indicate objective and subjective nights. 
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HvC02 only under LL, while the loss-of-function allele 
of SbPRR37 abolished the second evening peak of SbCO 
under LD and LL, as compared to a bimodal expression 
of SbCO in the Sorghum wild type. Taken together, our 
results suggest that PRR37 is regulated differently in long 
and short day grasses and mutations in the gene have 
differential effects on downstream genes, photoperiod 
sensitivity and flowering. 

Interestingly, we identified strong correlations in the 
expression patterns of HvC01/HvC02 and Vrn-Hl 
under LD and LL. Correlation of expression patterns 
may indicate that Vrn-Hl expression is affected by Ppd- 
Hl through HvC01/HvC02. Expression of Vrn-Hl in 
Scarlett was reduced, but cycled in LD, and was dam- 
pened with a trend towards arrhythmicity in LL 
(Figures 4g, h). It has been reported that Vrn-Hl is pri- 
marily controlled by vernalization, while photoperiod 
affects Vrn-Hl expression only indirectly through con- 
trolling Vrn-H2 expression [32] . Our results indicate that 
Vrn-Hl expression is circadian controlled and directly or 
indirectly affected by the mutation in Ppd-Hl. In Arabi- 
dopsis, the Vrn-Hl orthologs API and FUL have not 
been reported to cycle under constant conditions [69]. In 
barley, HvVRT2 and HvBMl, which act as repressors of 
flowering in barley [31], were also under circadian con- 
trol. Both genes showed higher expression levels in the 
late flowering genotype Scarlett, but diurnal and circa- 
dian cycling did not differ significantly between the two 
genotypes. Lower expression of HvVRT2 and HvBMl in 
the early flowering S42-IL107 {Ppd-Hl) than in the late 
flowering line Scarlett {Ppd-Hl) may thus reflect differ- 
ences in development rather than direct effects of the 
mutation in Ppd-Hl, 

Conclusion 

It was first shown in wheat that transcript accumulation 
of the Cab-1 gene, encoding the light-harvesting chloro- 
phyll a/b binding protein, was under circadian control 
[54] . However, the more tractable model plant Arabidopsis 
served then to advance the mechanistic understanding of 
how the circadian clock keeps time in plants. Recent tech- 
nical advances enable us now to return to the crops and to 
study inter- and intraspecific modulation of circadian 
time-keeping and its effects on adaptation and fitness in 
crops. The present study has demonstrated a high degree 
of conservation of the circadian clock genes in barley as 
compared to Arabidopsis and rice. However, differences in 
gene numbers, diversity, and in the function of Ppd-Hl 
suggest evolutionary modification of clock related compo- 
nents. Our results indicate that Ppd-Hl in barley 
has evolved novel functions in the control of flowering 
time, which are independent of its function in the circa- 
dian clock. Direct or indirect regulation of Vrn-Hl by the 
circadian clock and Ppd-Hl suggest modification of the 



photoperiod response pathway in barley as compared to 
Arabidopsis where the orthologs API and FUL are not 
known to be under circadian control. In addition, our 
results suggested differences in the regulation of PRR37 
{Ppd-Hl) in the long-day grass barley as compared to the 
short-day grass Sorghum. Structural and functional 
characterization of the barley circadian clock will set the 
basis for future studies of the adaptive significance of the 
clock in Triticeae species. 
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